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Abstract. We report detection of moderate to high-mass star formation in an infrared dark cloud (Gil. 11-0. 12) 
where we discovered class II methanol and water maser emissions at 6.7 GHz and 22.2 GHz, respectively. We also 
observed the object in ammonia inversion transitions. Strong emission from the (3,3) line indicates a hot ( 60 K) 
compact component associated with the maser emission. The line width of the hot component (4 kms -1 ), as well 
as the methanol maser detection, are indicative of high mass star formation. To further constrain the physical 
parameters of the source, we derived the spectral energy distribution (SED) of the dust continuum by analysing 
data from the 2MASS survey, HIRAS, MSX, the Spitzer Space Telescope, and interferometric 3mm observations. 
The SED was modelled in a radiative transfer program: a) the stellar luminosity equals ~1200 L© corresponding 
to a ZAMS star of 8 Mq; b) the bulk of the envelope has a temperature of 19 K; c) the mass of the remnant 
protostellar cloud in an area 8 x 10 17 cm or 15" across amounts to 500 Mq, if assuming standard dust of the 
diffuse medium, and to about 60 Mg, should the grains be fluffy and have ice mantles; d) the corresponding 
visual extinction towards the star, Ayr, is a few hundred magnitudes. The near IR data can be explained by 
scattering from tenuous material above a hypothetical disk. The class II methanol maser lines are spread out in 
velocity over 11 km/s. To explain the kinematics of the masing spots, we propose that they are located in a Kepler 
disk at a distance of about 250 AU. The dust temperatures there are around 150 K, high enough to evaporate 
methanol-contaiiimg ice mantles. 
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1. Introduction 

Infrared dark clouds (IRDCs) are cold, dense molec- 
ular clouds seen silhouetted against the bright diffuse 
mid-infrared emission of the Galactic plane. They were 
discovered during mid-infrared imaging surveys with 
the Infrared Space Observatory (ISO) ijPerault et alJ 
Il996tl and the M id-course Space Experiment (MSX 1 ) 
ijEean et alll 9981. The mm/submm observations revealed 
that typical IRDCs have gas densities of n > 10 6 cm~ 3 , 
temperatures of T < 20 K, and sizes of 1-10 pc 
ijCarev et alJll99S|h To determine the mass distribution 
toward the IRDCs, Carey et al. (2000) obtained images of 
the 450 and 850 /im emission using SCUBA. All observed 
IRDCs contain 1-4 bright sub-millimeter emission peaks 
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(> 1 Jy/beam at 850 fim) surrounded by an envelope of 
emission, which matches the morphology of the IRDC in 
mid-infrared extinction. 

ICarev et all ( 2000l) suggest that IRDCs might be sites 
of early phases of high-mass star formation, however stud- 
ies of their star formation content are still rare. There 
are, indeed two cases where detailed studies were done 
towards the IRDCs, although the original object identifi- 
cation was not based on the extinction in the MIR. Based 
on the submm dust continuum and line observations, 
ISandelll l|2000|) conclude that NGC6334 I ( N) is a high- 
mass Class object. Similarly. iGarav et al.1 1|2002|) suggest 
that IRAS16272 is a dense massive core in a very early 
evolutionary stage, distinguished by being luminous with- 
out any associated Ultracompac t Hii regions (UCHi iR). 
Toward the IRDC G79.3+0.3 Pl. lRedman et al.ll|2003|) re- 
ported associa t ed low and intermediate mass YSOs. Also, 
iTevssier et alJ l)2002l) found young embedded stars asso- 
ciated with the IRDCs they studied, but do not give fur- 
ther details on th em. Observations of IRDC Gil. 11-0. 12 
are discussed by Ijohnstone et alJ (2003) and, although 
they noticed signs of star formation associated with their 
SCUBA source PI, they did not find evidence of massive 
star formation. Identifying the earliest stages in the for- 
mation of massive stars is currently a crucial step in our 
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understanding of massive star formation, so we started a 
project to search for methanol masers toward cold massive 
cores. 

Class II methanol (CH 3 OH) masers (CIIMMs) are 
almost always found in high-mass star-forming regions 
(HMSFRs). It turned out th at not all of the CIIMM 
are associated with UCHiiR ijWalsh et al.lll998(l . Based 
on high spatial-resolution radio continuum and 6.67 GHz 
methanol spectral-lin e data in 364 s ources towards 
IRAS-selected regions IWalsh et all l|l998l) show that the 
methanol maser is most likely present before an observable 
UCHiiR is formed around a massive star and is quickly 
destroyed as the UCHi i R evo lves. However, this is dis- 
puted by IPhillips et alJ (|l998^) who argue that the non 
detection of radio continuum emission associated with 
methanol masers is because of the association of the 
maser sites with low/intermediate mass stars which do not 
ionise their immediate surroundings. To date no CIIMM 
has ever been detected toward a low-mas s SFR. with 
the p ossible exception of NGC 2024-FIR4 l|Minier et alJ 
2003), which, by the way, is the only known CIIMM in 
the whole O r ion G iant Molecular Cloud (GMC) complex. 
IWalsh et al.l l|200l[) in another study of 31 methanol maser 
sites and 19 UCHiiR in Mid-infrared (10.5 and 20 /xm) 
find that the observations are consistent wi th the maser 
emissi on being powered by the MIR source. I Walsh et alJ 



(2001) also find stron g evidence of in f all/ou tflow in differ- 
ent molecular tracer. iBeuther et alJ l|2002l) find a mean 



separation of ~ 0.19 pc between the CH3OH and cm 
emission for there sample of 29 massive star forming re- 
gions. Subsequently , in a continuum mapping reported by 
IWalsh et all l(2003ti toward 71 6.7 GHz maser fields, they 
find that most methanol masers are within 10" of a sub- 
mm peak. Furthermore, a recent 1.2 mm continuum survey 
in massive star forming regions with methanol m aser sites 
and /o r radio continuum emission, reported by IHill et alJ 
( 2005) shows that these tracers of massive star formation 
are associated with heavily embedded protostars. These 
arguments taken together with the non detection of radio 
continuum make a convincing case of CIIMMs signposting 
the earliest stages of high mass star formation. 

Here we report the discovery of CIIMM emission to- 
ward the IRDC Gil. 11— 0.12 and discuss the properties 
of the protostar heavily embedded in the IRDC detected 
with SPITZER in the MIR. In ® we describe our ob- 
servation in various molecular tracers with the Effclsberg 
100-m telescope, Australian Telescope Compact Array, 
Very Large Array and the Berkeley-Illinois-Maryland- 
Association interferometer. In Sjlllwe present the evidences 
for active star formation in the IRDC and model the am- 
monia emission. The discussion on the Spectral Energy 
Distribution (SED) of the embedded source and modelling 
in spherical symmetry follows through in [JH The observed 
properties of the CIIMM are entirely comparable with 
the properties of numerous other such sources detected 
in high-mass star forming regions. Therefore, its presence 
makes Gil. 11— 0.12 , an IRDC in which active high-mass 
star formation has been unequivocally established. 



2. Observations 

2.1. CH 3 OH maser 

We observed the class II methanol maser at 6.7 GHz 
(5i — 60 A + transition) toward the submm peak IRDC 
Gil. 11 PI (a J200 o = 18 : 10 : 28.402, £,2000 = -19 : 22 : 
29.00) with the Effelsberg 100-m telescope 2 on 15 June 
2003. The frontend was the facility 5 cm receiver tuned to 
6668.518 MHz. We used the 8192 channel auto-correlator 
with 2 subunits of 10 MHz bandwidth to be able to reach 
a spectral resolution of 0.25 kms -1 . The beam width of 
122" was determined from drift scans over compact con- 
tinuum sources, which also served to check the pointing. 
The observations were repeated with the 100-m telescope 
on 5 August 2004. 

On 11 January 2004, we followed up the observa- 
tions of the 6.7 GHz Class II CH3OH maser with the 
Australian Telescope Compact Array (ATCA) in the 
FULL_4_1024-128 correlator configuration, employing the 
snapshot mode with a series of 6 short integration cuts 
of 5 minutes. Since ATCA can observe two frequencies 
simultaneously, we chose the 8.64 GHz continuum band 
covering a bandwidth of 128 MHz frequencies. The 1 a 
sensitivity achieved is ~ 0.2 mJy/beam. Phase calibration 
was done using B1817-254. The primary beam of ATCA 
at this frequency is ~ 8 and the effective HPBW of 
the synthesized beam equals 2.03" x 6.24". The spectral 
resolution was 0.2 km/s, and the 1 a sensitivity 0.5 Jy. 
The ATCA spectrum is essentially identical to the 2003 
Effelsberg spectrum. 

2.2. Ammonia observation 

The maser detection prompted us to search for the NH3 
(3,3) and (4,4) transitions with the Effelsberg 100-m tele- 
scope in April 2004. With the AK 8192 backend, we were 
able to observe the (1,1), (2,2), (3,3), and (4,4) transitions 
in both polarisations simultaneously using the K-band re- 
ceiver. With 8 subunits of 10 MHz bandwidth, the result- 
ing spectral resolution was ~ 0.2 kms -1 after smoothing 
the data to improve the signal-to-noise ratio. The beam 
at the NH3 frequencies was 40". The observations were 
performed in the frequency switching mode. Pointing was 
checked at hourly intervals by continuum scans on G10.62. 
We estimate, the pointing to be accurate within 6", with 
the pointing scans used for the absolute calibration. 

2.3. VLA H2O maser observation 

The 22.2 GHz H2O maser observations were done with 
the Very Large Array (VLA) on 24 August 2004 in its D 
configuration and in 2 polarisations with a spectral reso- 
lution of 0.33 km/s, a synthesized HPBW of ~ 0.4", and 
a primary beam of ~ 120". The standard interferometer 
mode was used with a total integration time of 20 minutes 



2 Based on observations with the 100-m telescope of the 
MPIfR (Max-Planck-Institut fr Radioastronomie) at Effelsberg 
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on source, split into sessions. Phase calibration was done 
using J1733-130, and the flux calibrator was 3C286. 

2.4. BIMA observation 

The 3mm observations were carried out from September 
1999 to July 2000 with the Berkeley-Illinois-Maryland- 
Association interferometer (BIMA) in its B, C, and D 
configuration sampling spatial structures from 3" to 60". 
Continuum images were obtained with a combination of 
86 and 93 GHz tracks centred with a beam size of 8.3" x 
3.9" reaching an rms of 1.6 mJy The total bandwidth of 
both sidebands was ~ 900 MHz. 



3. Results 

3.1. CH 3 OH maser 

The solid line in FigQ] displays the 6.7 GHz methanol 
maser spectrum toward G11.11P1 detected in August 
2004 with the Effelsberg 100-m telescope, while the filled 
region shows the June 2003 data. As the Effelsberg 100- 
m absolute calibration has an uncertainty of ~ 20%, the 
observations taken with two different instruments and at 
different times have been normalised with respect to one 
of the features, so that absolute calibration is no longer an 
issue. The maser features cover a velocity range of about 
llkm/s, with the peaks at 32km/s in June 2003 and at 
24km/s in August 2004. Thus there is clear evidence of 
a variation in the relative line intensities of about 10 %. 
From mapping the source, we find that the masing spots 
originate at the position PI from an area that is point- 
like with respect to the Effelsberg beam. There is no evi- 
dence of another maser source. The follow-up observations 
obtained six months later with ATCA reveal no obvious 
changes in the spectrum and also indicate that the maser 
originates from a single core. The peak flux of the bright- 
est maser feature from the ATCA observations is 22 Jy. 

The position of the integrated maser emission is plot- 
ted in Fig 2 There, the left panel shows the MSX map of 
the filame ntary dark cloud ove rlaid on the SCUBA submm 
map from ICarev et al.1 l|2000h . and the right panel is a 
blow-up of the SCUBA image around PI, together with 
contours of 2MASS K s 3 band and the 3mm continuum 
emission. To estimate the positional accuracy of the maser 
observations, we split the phase calibrator observations 
into two sets and cross-calibrate the second set with a self- 
calibrated first set, which results in a position accuracy of 
(0.05",-0.15"). This is a lower limit, but since our source is 
less than 10° away from the calibrator, baseline errors are 
not likely to add to the error. To cross-check this, we re- 
duced other observations towards an already known maser 
source and found a positional accuracy of (0.105", -0.75"). 

3 The atlas image obtained is part of the Two Micron 
All Sky Survey (2MASS), a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the 
NASA and the NSF 
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Fig. 1. Lower panel: Effelsberg 100-m Spectra of the 
CH3OH emission at 6.7 GHz. The filled spectrum cor- 
responds to June 2003 and the solid lines to August 2004 
observations. The blue shifted component is normalised to 
unity in both cases. Upper panel: 22.2 GHz water maser 
spectra obtained with the VLA. 
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Fig. 3. Velocity structure of the CH3OH maser emission. 
Each data point gives the position of the model fits to the 
individual velocity channels. The corresponding Vlsr m 
kms _1 is also shown. 



ATCA is an east west array, so the beamwidth in decli- 
nation is greater by a factor cosec(DEC); and for those 
sources with declination north of -24 deg, complete u-v 
coverage is unobtainable. 

When the sub-spots are analysed in position-velocity 
plots, they fall nicely along an arc, as predicted for a disk 
l|Minier et al.ll998llNorris et alJll998|) . There is also a lin- 
ear velocity trend. This is shown in Fig|3] where a point 
source model was fitted to individual velocity channels. 
The o ne sigma position uncertainty is proportional to the 
SNR (jReid et alJll988|) : hence, with a SNR 15-30 for the 
different velocity features, we reach a sub-arcsecond ac- 
curacy. The relatively weak peak at 34 km/s was not in- 
cluded due to its large positional uncertainty. 
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Fig. 2. Left: The 8 jum image of Gil. 11 with SCUBA 850 /mi (Carey et al. 2000) overlay Right: BIMA 3mm continuum 
image towards PI in contours, with contour level~1.6 mJy (-3,3,5,7) and 2MASS K s band in dashed contours on the 
850 (im image. The star denotes the H2O maser and the filled circle the CH3OH maser position. 



3.2. 3mm continuum emission 

The BIMA image shows a compact, slightly resolved 3mm 
continuum source that coincides with the SCUBA sub-mm 
peak shown in Fig|2 A 2-D Gaussian fit to the emission 
yields a source size of0.25x0.12pc with a PA 3° . The peak 
flux of ~ 12 mJy corresponds to a brightness temperature 
of ~ 28 mK. Our cm continuum observations with ATCA 
at 8.4 GHz, along with searches in Galaxy-wide cm contin- 
uum surveys, reveal no detectable free-free emission. In or- 
der to check whether one expects to find free- free emission 
at 3mm, we extrapolated the ATCA upper limit to that at 
3mm, assuming optically thick emission. We estimated a 
mm flux of 25.3 mJy relative to the ATCA upper limit of 
~ 0.2 mJy. Thus the contribution from free- free emission 
cannot be ruled out. For a typical value of 10000 K for the 
ionisation temperature from an optically thick ionised gas 
cloud, we get an emission measure EM > 4 x 10 10 pc cm" 3 
at 90 GHz , a very high value that is more suit ably as- 
signed to Hypercompact H11 regions l|Kurtzl 2002'). Such a 
region with typical sizes < 0.05 pc would have been point- 
like with respect to the BIMA beam, while we found that 
the source is partially resolved. Thus, the 3mm emission 
is most likely due to dust. 



3.3. H2O maser 

The position of the H2O maser appears offset from the 
CH3OH maser position, as well as from the BIMA 3mm 
continuum peak, within the respective positional errors. 
This might indicate that, although both CH3OH and H2O 



masers occur in warm and dense environments, they may 
not sp atially coincide due to different excitation mecha- 
nisms feeuther et al.ll2002j) . 



Water masers can be explained by collisional pump- 
ing with H? molecu les in shocks associated with outflows 
l|Elitzur et alJ lT989). Alternativel y, the pumping can oc - 
cur by accretion shocks in disks ijGarav fc LizanolH9 991. 
The maser is weak (~ 0.3 Jy) with a feature offset of 
~ 7 kms _1 from the systemic velocity as shown in the 
upper panel of Fig^ This agrees well with the picture 
of water masers being variable and having large velocity 
spread. 



The maser spot is not spatially resolved over its ve- 
locity structure. We performed a model fit to the indi- 
vidual velocity channels over the single feature. Unlike 
the CH3OH masers, the H2O maser spots do not show 
any systematic gradient. Since the velocity structure is 
also found to be spread on a much higher angular scale 
than the CH3OH maser, it may be tracing the outflow. 
Higher angular resolution is needed to reveal structures 
that could be either a shock front of the high-velocity 
outflow very close to the central object or part of the 
accretion disk itself. We also serendipitously discovered 
another water maser that was offset by more than 1'. The 
feature is brighter and situated at the systemic velocity of 
the source. The absolute position of the maser peak is at 
aj 2 ooo=18:10:33.58,5j 20 oo=-19:22:50.3. 
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3.4. Ammonia towards G11.11P1 

Figure 01 shows the ammonia lines observed with the 100- 
m telescope toward G11.11P1, and the line parameters are 
given in Tabled F° r the (4,4) line that was not detected, 
the ler r.m.s noise in the spectrum is given instead. 

Table 1. Effelsberg 100-m line parameters 



Transition 


VLSR 


Tmb 


FWHM 




(kms -1 ) 


(K) 


(kms -1 ) 


NH 3 1-1 


29.78(0.01) 


4.1 (0.4) 


1.29(0.02) 


NH 3 2-2 


29.78(0.02) 


1.58(0.07) 


1.81 (0.06) 


NH 3 3-3 


29.77(0.07) 


0.64(0.05) 


3.22 (0.21) 


NH 3 4-4 




(0.065) 





LnH 3 (4,4) z 


^NH 3 (3,3) J 


- NH,(2,2) 

i J 


{ \ 


r nh 3 ( 1,1) 





40 : 

Velocity (km/s) 



Fig. 4. top to bottom: Spectra of the NH 3 (J,K) = 
(1, 1), (2, 2), (3, 3), and (4,4) transitions. All spectra were 
taken towards the submm peak position given in the text. 
The fit discussed in *13.4l is overlaid in red. 



nication). The program then yields the rotation tempera- 
ture, column density and source size for a given telescope 
size assuming a Gaussian structure, and produces a syn- 
thetic spectrum in agreement with observations. Table 2 
lists the fit results. Thus the broader component appears 
to be much hotter (60 K) and is very compact (3.3"). 
Actually, the fit is not unique since there is degeneracy 
between column density and the source size. But a large 
source size would mean a still higher luminosity, which 
would put the source in the spectral type 09. It is highly 
unlikely that an HII region around such a star goes unde- 
tected. On the other hand, a much smaller source size with 
very high column density would mean extremely dense gas 
with critical density high enough to excite the (4,4) level, 
which we do not detect. 

Additionally, we fi t the ammonia spectra with a spher- 
ical symmetric model ijWvrowski et al.ll2003ft using both a 
core component with a temperature power law, consistent 
with internal heating from a 1200 L© embedded source, 
and a more extended, cooler (10 K) lower density clump. 
For the core component, the luminosity, and the outer and 
inner radii were taken from the dust continuum modelling, 
which will be presented in i)4.3.1l Assuming an ammonia 
abundance of 10 -8 , the density in the core is 7x 10 5 cm~ 3 , 
which is consistent with the modelling in H4.3.1I and is a 
factor 10 lower in the cooler clump with a size of 40". For 
the core component, a line width gradient was modelled 
with 4 km/s in the inner 2 arcsec and then decreasing 
to the clump component with 1.2 km/s. The increasing 
line widths/turbulence towards the centre of the core is 
another indication of an embedded, massive young stellar 
object that is stirring up its environment, as well as the 
remaining line wings of the (3,3) line, which indicates an 
additional outflow component. 



Table 2. Model fit results 



Component 


Trot 


size 


-/Vnh 3 


FWHM 




(K) 


(") 


(10 15 cm" 2 ) 


(kms" 1 ) 


NH 3 cold 


15.4 


23 


1.8 


1.34 


NH 3 hot 


60 


3.3 


13.7 


4 



The line widths from Gaussian fits when taking the 
hyperfine satellites for the (1,1) and (2,2) spectra into ac- 
count increases with the excitation of the lines. The (3,3) 
transition even shows, a broad wing besides the line core. 
From a Boltzmann plot of the ammonia emission, we con- 
clude that the ammonia emission cannot be explained with 
a single temperature. From just the NH3 (1,1) and (2,2) 
data, we got a rotational temperature of ~ 14 K. The ro- 
tational temperature derived from the (2,2) and (3,3) lines 
was estimated to be ~ 30 K. 

This evidence of two temperature components 
prompted us to fit the spectrum with two components si- 
multaneously using XCLASS (Schilke, P., private commu- 



3.5. Archival Data 

We collected the data from 4 galactic sky survey archives, 
namely 2MASS, Spitzer, MSX, and IRAS. The 2MASS 
point source catalog contains the JHK fluxes for 3 point 
sources all within the 14" of the SCUBA beam. The 
Spitzer MIR data collected from the Spitzer data archive 
is part of the GLIMPSE Legacy survey in 4 of the Infrared 
Array Camera (IRAC) bands from 3.6 — 8/im. We found 
3 point sources in the 3.6/^m band, while in the higher 
wavelength bands we found a single source that were co- 
incident with the MSX 8/mi point source. The FIR emis- 
sion obtained from the high resolution IRAS data is not 



6 



Pillai et al.: high-mass star formation in Gil. 11 



point like and might suffer from background emission. But 
IRAS point sources generally do not appear circular, pri- 
marily because the detector aperture were rectangular. We 
adopted an flux uncertainty of 20% in all cases as given 
by the respective Space centres. 

4. Discussion 

4.1. CH3OH maser disk scenario 

Figure ^ contains 3 major peaks that is very reminiscent 
of a classical maser triplet, which is usually ascribed to 
a maser amplifi cation from a disk in Kep lerian rotation 
viewed edge on llElmeereen fc Morrislll979(l. 

Ac cording to ICesaronl ~ i^)90l ) and IPonomarev et al.1 
l)l994|) . the intensity ratio of the central to the red and blue 
peaks is determined by the width of the ring in which the 
conditions for line amplification are favourable, and the 
central peak vanishes when the ring becomes very small. 
The red and blue components of our spectrum are not 
equidistant to the central component that marks the sys- 
temic velocity (see also Table[3]of NH3 velocities). Such a 
distortion may be due to a tilt of the disk, in which case 
the physical parameters (density, temperature, and veloc- 
ity) along the line of sight depend not only on the radial 
distance, but also on the height above the disk. 

In principle, equidistant triplets can also arise from 
an expanding shell. We favour the disk scenario for sev- 
eral reasons, a) Variations in the line ratios, as argued by 
ICesaroml l|l990l) . are readily explained for disks by radia- 
tive interaction between masing spots, whereas such an 
interaction would not be possible for an expanding shell 
because of opposite velocities of the masing. b) All 3 ma- 
jor peaks contain sub-peaks, which we interpret as arising 
from irregularities associated with the tilt or warp of the 
disk, c) The near-infrared emission and the Spitzer data 
discussed in §4.1.2 also seem to require a disk structure. 

4.2. Mass & luminosity estimates 

The kinematic distance to Gil. 11 derived as 3.6 kpc us- 
ing the rotation curve of Clemens (1985) and the IAU 
standard values of Rq — 8.5 kpc and Vn = 220 k ms^by 
Carey et al. (1998) is 3.6 kpc. lCarev et alJ l|200(Jl quote a 
gas mass of 67 M Q from their dust continuum observations 
at 850 micron. The virial masses estimated from ammo- 
nia observations range from 150 — 240 M Q for the hot and 
the cold components within the beam (40") of Effelsberg. 
Similarly, by assuming the fractional abundance of NH3 
(relative to H2) to be on the order of 10~ 8 , the NH3 col- 
umn density translates to masses of 200 — 1600M Q . 

The CIIMM need elevated temperatures for two rea- 
sons: first, high methanol abundances are required to pro- 
duce an adequate maser gain,, and second, high tempera- 
tures are needed to produce pump photons. To effectively 
evaporate methanol off icy grain mantles, temperatures 
in excess of 100 K are needed (Sandford & Allamandola 
1993). Recent models of CIIMM pumping require gas tem- 



peratures ~ 100 K, somewhat higher dust temperatures, 
and densities of about 10 7 cm~ 3 . Therefore, the compact 
CH3OH maser emission we observe must arise from a hot, 
dense source. We may adopt the spread of the maser com- 
ponents, 0.25" or 900 AU, as the size of that source, which 
is comparable to the sizes of numerous other CIIMM re- 
gions. The Stefan-Boltzmann law delivers with a luminos- 
ity of 830 Lq for a source of the corresponding area and 
T = 100 K . To fit all observed ammonia lines simulta- 
neously, we presented a model in H3.4I that assumes an 
optically thin spherical dust cloud illuminated by a cen- 
tral point source, which then results in a radial tempera- 
ture distribution of T cx r" 4 (Wilner et al. 1995). With 
a luminosity of 1200 Lq, we could explain the observed 
warm ammonia core embedded within a larger cold clump. 
In the model, the temperature reaches a temperature of 
240 K at the inner radius, hence large enough to excite 
the methanol maser observed. 

4.3. SED 

4.3.1. MIR and longer wavelengths 

The SED of Gil. 11 PI shown in FigEl contains fluxes due 
to dust emission at (i) 3mm (BIMA, 8" x 4" beam) , (ii) at 
450 and 850/im (JCMT), (iii) IRAS data, (iv )MSX upper 
limits at 12 and 21 /mi, (v) an MSX detection at 8.2 /im, 
(vi) four detections with the Spitzer Space Telescope (at 
3.6, 4.5, 5.8 and 8.0/im), and (vii) three NIR data points 
from 2MASS. 

Although for massive young stars, a time sequence 
from class to class 3 has not been established, a high- 
mass-star equivalent to a class object should also emit 
more th an half a percent of its energy at submm wave- 
lengths l)Andre et al.lll993t) . For Gil. 11, the luminosity at 
A > 300/im amounts to almost 3% of the total, so in this 
respect it could be considered an object equivalent to a 
Class source. Similar arguments of a Class equivalent 
object in hi gh-mass-star - formin g region have also been 
proposed bv lMotte et alJ 1)2003|) in the submm protoclus- 
ters in W43 . They derive a submillimeter to bolometric 
luminosity ratio of 1.5-3% in W43-MM1 and MM2. 

To estimate how the dust temperature and density 
vary within the envelope, we tried to reproduce the con- 
tinuum emission in a self-consistent way by calculating the 
radiative transfer in spherical symmetry. However, when 
including the NIR fluxes, as well as the recently published 
data from the Spitzer archive, the assumption of strict 
radial symmetry had to be abandoned. We sketch below 
a possible, or even likely, structure of the source; but to 
confirmat its correctness we need information with higher 
spatial resolution, a difficult requirement for such a distant 
source. The model therefore remains rather qualitative. 

Nevertheless, the observations yield a number of im- 
portant constraints, a) The total luminosity is determined 
by the IRAS fluxes. While the 100/im flux may contain 
a fair contribution by diffuse emission (from cloud heat- 
ing by the interstellar radiation field or scattered stellar 
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sources), the 60/im flux indicating warm dust comes from 
a compact source with a luminosity of ^10 3 L©. b) As the 
source is only marginally resolved at 3 mm, its outer cloud 
radius r out has to be smaller than 4 x 10 17 cm. c) The 
submm/mm fluxes lead to a total mass of several hun- 
dred M when adopting standard interstellar dust with 
an absorption coefficient a t 1 mm of ~4 x 10 -3 cm 2 per 
gram of interstellar matter l)Krugell2003j) . However, as this 
emission is optically thin, a substantial part may not be 
directly linked to the star, but may come from the back- 
ground within the IRDC. d) The visual extinction to the 
star depends on how the dust is distributed. It is about 
400 mag if all the mass is in a sphere with constant density 
around the star, and even higher for a centrally condensed 
envelope, and less than 400 mag in case of an additional 
background cloud. For comparison, we derive a column 
density Ay > 60 mag from the peak of the ammonia emis- 
sion in a 40" beam. This column density is probably two 
or three times larger, depending on the compactness of the 
source. We then assumed conversio n factors [NHal / [H o] = 
10" 8 and A H = 2 x 10 21 Ay cm" 2 llMotte et nl.lll!)!),4. 

The solid line in FigO which is compatible with the 
observations between 10/im and 3mm, comes from a spher- 
ical radiative transfer model (see Appendix A) with a cen- 
tral star of L = 1200 L Q and 15000 K surface temperature 
that is surrounded by a envelope whose inner and outer 
radii are 10 15 and 3.5 x 10 17 cm, respectively. The exact 
number for T e g is irrelevant, because the stellar radiation 
gets immediately absorbed. The dust density in the enve- 
lope is constant (pdust = 3 x 10 -20 g cm -3 ) implying an 
optical extinction Ay ~ 400 mag. It may, at most, be two 
times less, otherwise the model spectrum would not com- 
ply with the 21^,m upper limit. If it is less than 400 mag, 
one has to invoke the above mentioned background cloud 
to explain the sub- mm fluxes. These model results are also 
consistent with our fits to the ammonia lines described in 

an 

It is important to note that the model spectrum for 
the FIR does not rise again at shorter wavelengths and 
cannot account for the 8.2/xm MSX and the Spitzer fluxes. 
Here an independent source is required. We propose that 
besides the star, there is a second blackbody emitter close 
to it. Its temperature is approximately T — 550 K, its 
total luminosity 150 L©, and it suffers 150 mag of visual 
extinction. This is below the minimum estimate for Ay 
of the envelope, but the latter may, of course, be patchy. 
It is natural to associate the second blackbody with an 
accretion disk of about 12 AU radius (outer). Considering 
the paucity of the data, more refined calculations seem, at 
present, unwarranted. 

The dashed line that fits the long wavelength data in 
FiglSJis a modified Planck curve of the form v 2 B v (T) with 
T = 19 K, possibly the emission of the background cloud 
with a total mass of 350 M© . 




1 10 1 10 2 10 3 

wavelength [/U.m] 



Fig. 5. The Spectral Energy Distribution of Gil. 11-0.12 
PI. Squares with errorbars indicate observations, where 
two almost coincide at 8 and 8.2/im. The upper limits from 
MSX corresponding to bands C (12.13/im ), D (14.65/xm), 
and E (21.3//m) are shown. The solid curve represents the 
semi-spherical model described in the text. The solid line 
at FIR wavelengths corresponds to the SED of a central 
star of L = 1200 L© and to the 15000 K surface temper- 
ature that is surrounded by an envelope. Possible second 
blackbody emitter is shown as solid curve at shorter wave- 
lengths. The dashed curve represents a grey body fit of 19 
K. The NIR fluxes is explained by scattering (see H4.3.2|) 

4.3.2. 2MASS sources 

As shown in Fig[3 the CH3OH peak is offset from the 
2MASS source by ~ 2". It is highly unlikely that the 
ATCA observations have errors on the order of 2" (cf §3.1). 
Though the rms uncertainty in the case of 2MASS data is 
on the order of 0.4", we cannot completely rule out that 
the embedded object triggering maser emission is not as- 
sociated with the 2MASS source. Because of the high col- 
umn density towards the centre of the core, it could be 
that we are seeing light from the embedded object scat- 
tered out from the edges of the NIR core. We now discuss 
such a situation in detail. 

The 2MASS J and H bands reveal 3 faint point sources 
detached from the sub-mm peak, while at K, where the 
resolution is inferior, there is only one single source. These 
2MASS points lie far above any possible fits to a spher- 
ical cloud model. We first investigated whether they can 
come from an unrelated foreground star of lower bolomet- 
ric luminosity but suffering less obscuration. Analysing 
the observed JHK fluxes in a colour-colour diagram, we 
found that they cannot be explained by reddening, so we 
discarded the possibility of a foreground star. 

Instead, we propose that we are dealing with three 
knots of scattered light that escapes from the star into an 
optically thin cone above a hypothetical disk. We estimate 
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the flux scattered by the cone at frequency v from the 
formula 



S,j — 



4ttD 2 ' 



(1) 



where D is the source distance and L v the stellar spec- 
tral luminosity, so L — J L v dv. If f2 C one denotes the solid 
angle of the cone and r^ ca its (small) scattering optical 
depth, then /„ = r;f a fi cone /47r. Because of the basically 
unknown geometry, our estimates are very rough. But as- 
suming the star to be a blackbody, only about one percent 
of the light has to be scattered (the exact number depends 
on the stellar temperature). 

As regards the wavelength dependence, scattering of 
pure stellar light would suggest decreasing fluxes at longer 
(K) wavelengths. However, foreground extinction and, 
more importantly, hot dust from an accretion disk (T ~ 
1000 K) could easily explain the observed spectral shape at 
near IR wavelengths. In view of the geometrical require- 
ments posed by the near IR fluxes, we propose a more 
realistic configuration. In such a model, the protostar is 
part of a massive (several 100 M©) clumpy cloud com- 
plex. Such a patchy dust distribution is shown in Fig|B| 
As seen in Fig El the emission at NIR is concentrated at 
only one edge of the BIMA 3mm emission. In the configu- 
ration envisaged in FigEI the stellar radiation is scattered 
by the dust clump above the disk at NIR wavelengths to- 
ward the observer. However, the stellar photons scattered 
from regions below the disk suffer severe extinction from 
the disk and fail to reach the observer at NIR wavelengths, 
resulting in a uni-polar nebula in NIR. 

We therefore regard the JHK knots as further indirect 
evidence of a stellar disk. 



5. Conclusion 

Our main result is the detection of what is most likely a 
high-mass young stellar object within the IRDC G 11.11. 
If the luminosity we derived is an overestimation that cor- 
responds to an intermediate star, this would be the first 
reported case of CIIMM detection towards an IM mass- 
star-forming region. This evidence comes from the follow- 
ing observations and arguments: 

— CH3OH and H2O masers appear during the earliest 
phases of massive star formation, so their detection in 
Gil. 11 is evidence of an embedded object. 

— The systematic velocity gradient and the linear feature 
in the P-V diagram of the CH3OH lines is a signature 
of a disk encompassing a protostar. 

— The higher excited NH3 (3,3) thermal line points to- 
wards a compact and hot central component. The ther- 
mal CH3OH emission (Leurini et al. in preparation) 
suggests an outflow, in agreement with the claim by 
I Johnstone et alJ ll2003l). based o n H2CO obervations. 
However. I.Tohnstone et all l)2003h do not find any evi- 
dence of high-mass star formation. 



— The strong FIR fluxes (> 60 Jy) derived from 
analysing the high resolution IRAS data imply a pow- 
erful central source (> 1000 L©). 

— The radiative transfer additionally requires blackbody 
emission of about 550 K from a separate, less embed- 
ded source, probably from a disk, in order to repro- 
duce the Spitzer data between 3.6 and 8.0^m. A non- 
spherical geometry is also indicated by the 2MASS 
NIR fluxes, which can only be explained by scatter- 
ing and not through reddening. 

— The high-extinction estimates also suggest an ex- 
tremely young object, probably the Class equivalent 
of a high-mass protostar. 

Appendix A: Model 

We iteratively solved the energy equation 



K v J v dv = J K, v B v {T A )dv 
together with two integral radiative transfer equations 



I+(r) = er T J+(0) + / S(x)e x dx 



r{t) = e -1 / S(x)e x dx 
Jo 

Here k v denotes the dust extinction coefficient, J„ the 
mean intensity, and B V (T&) the Planck function at dust 
temperature In the formulae for the radiative trans- 
fer, I + (t) is the intensity directed towards the observer 
at optical depth t, frequency is, and impact parameter 
p, while /~ refers to beams away from the observer. The 
source function S(t) includes a scattering term. The op- 
tical thickness r and geometrical length z are related 
through dr = pndz, where p is the density. Index v has 
been dropped for convenience of writing. 

There are two boundary conditions. One states that 
the standard interstellar radiation field (ISRF) impinges 
at the cloud surface (I~ = l\ , F ), the other reads 
I + — I~ = at r = (where also z = Q) at impact 
parameters p greater than the stellar radius i?* , and that 
I + — 1~ = F^v/ir at the stellar surface (F^y is the stellar 
flux at frequency v). 
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Fig. 6. Illustration of an intuitive model of the source geometry needed to explain the NIR emission. The protostar 
is part of a massive (several 100 M©) clumpy cloud complex. The optically thin outflow cone scatters the stellar light 
at NIR wavelengths toward the observer. The stellar photons scattered from the cone below the disk is extincted by 
the optically thick disk; consequently, we observe a uni-polar nebula at NIR. 



